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Analysis  of  navigation  systems  employing  independent  inertial  anti 
velocity  measurements ,  begun  in  Rel\  1,  is  continued.  Explicit  fcncolaa 
are  given  for  mlnivium  rms  error  in  computed  position  as  e.  function  (  f 
tiE9  of  flight.  Curves  based  on  these  formulas  are  presented,  showing 
the  results  for  a  nuiiiar  of  illustrative  cases. 
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LIST  OF  SYMBOLS 


Symbol 

lefinition 

Fage  First  Useu 

-'■■(t) 

vehicle  position  at  time  t 

1 

x(tjT) 

optimum  ccraputoris  estimate  of 
x(t)  based  on  all  the  dial 
readings  up  to  time  T 

1 

x(T) 

x(TjT) 

1 

T 

elapsed  time  since  beginning  of 
flight 

1 

Yl»  t2»  y3’  V 

parameters  describing  instrument 
error  statistics 

1,2,5,16 

°2 »  ’  uli  Pi 

6A(t}T) 

x(tjT)  -  x(t) 

2 

^x(T) 

x(T»T)  -  x(T) 

2 

defined  by  Eq,  (JI.l) 

7 

1J 

defined  by  Eq.  (in.l) 

15 

accelerometer  error 
autocorrelation  function 

l 

> 

CO 

w 

velocity  dial  error 
autocorrelation  function 

l 

12 

frequency  of  earth's  radius  pendulum 

3 

M- 

defined  by  Eq.  (iljli) 

It 
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1.  INTRODUCTION 

Tilda  report  is  an  extension  of  the  work  begun  in  Ref.  1.  This  is  not 
a  self-contained  report j  rather,  a  reading  of  Ref.  1  is  a  necessary  basi3 
for  the  understanding  of  the  following  rate  rial.  Repetition  of  expository 
material  and  results  obtained  in  lief.  1  will  be  kept  to  a  minimum. 

The  major  results  of  Ref.  1  were  obtained  for  two  canes ,  vizi 

Caso  1  0B(o, t)  » 

Case  c  0^(s,t)  ”  i"2 

where  00(s,t)  and  ^(s.t)  are  the  autocorrelation  functions  of  the  accelero¬ 
meter  dial  error  and  the  velocity  dial  error,  respectively..  Formulae  wera 
derived  for  the  optimum  method  of  position  computation  in  each  easoj  also, 
formulas  were  given  by  means  of  which  the  ras  error  in  computed  position  as 
a  function  of  time  of  flight  could  be  derived. 

In  the  following  pages ,  the  actual  formulas  for  rma  error  in  computed 
position  will  be  given j  also,  computations  based  on  these  formulae  will  bo 
given  for  a  number  of  illustrative  cases. 

As  in  Ref.  1,  let,  for  0  ■  t  ■  T, 
r(t)  »  true  position  at  time  t 

x(tjT)  *  optimum  computer's  estimate  of  x(t),  baaed  on  all  the 
dial  readings  up  to  tins  T. 

Denote  xlTjT)  simply  by  x(T) . 
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Also  let 


(IJ) 


M  t  jT,'  “  Jt(  t}T)  -  x(t) 

A 

■>(T)  »  ;^(T;T)  “  x(T)  -  x(T) 


and 


(1.2) 


c  ^ 

^(tjT) 


jx(tjT)  -  x(t)]‘ 
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the  mean  being  taken  over  an  ensemble  of  flights. 

Ae  in  Rof .  1,  tho  following  notation  in  used  for  description  of  the 
Hal  error  statists  cat 

Gn(t)  »  accelerometer  dial  error 

D 

t'y(t)  “  velocity  dial  error 

£_  “  independent  initial  position  dial  error 


£ 


Tl  “  Independent  initial  velocity  dial  error 

Ao 


The  en  a  amble  naans  of  all  these  errors  arc  zero)  also 


(1.3) 


t3 


and 


(I JO  -  £B(o)  O0(t) 


0y(a,t)  -  Gv(s)  £y(t) 

(All  averages  are  taken  over  an  ensemble  of  flights.) 
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II.  CASE  I:  ^B(a,t)  “ 


Lot.  3  “  lQ  t^  ' '  v  “  T  ho  n  equally  spaced  time  points  in 
the  interval  (0,1).  Let 


(II  .1) 


M 


1y(\,  tj) 


-1 


(matrix  inverse) 


Also  let 


n-1 

l>l.’nTL  Htj 


i,j“o 


(II.2) 


n-1 


^]n  "  0  HL  ‘id  8in  -°-tL  cos  Slij 


ij-o 


n-1 


LA22Jn  ‘  L _  5ij  COS  -Cl\  cos  ~Ltj 


i,J-o 


and 


An-  lia  jO 

a->oo  L  _ln 


(11*3)  1^2  ■  lim  j  A.  „ 

rl  \  _  I  “  ‘ 


n-»oo  -  in 


A22  “  :L1d  I A22  : 

n  -r>  o-  L  _  n 

The  quantities  A-^,  A-^,,  A^  are  functions  of  T, 
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ieft,rurf'e  1  gave  formulas  (Lqs.  111.36,  III. 17',  ITI.20  of  h~f.  1; 
from  which  could  be  der*  rod  the  formula  for  (L|(tjT).  Tho  result  can  bo 
expressed  as  follows : 


(II.!*) 


Let 


TF*3 

^*T3 


(ii. 5) 


and 


r 


Ln‘An  + 


\  *12  "  h.2 


^22  “  ^22 


rh 


A “  hi  *12  ”  *12 


(II. 6) 


Then,  for  0  ■  t‘ 

r  2  l 

C^UfT)  -  p  +  j 

/L'22 

L 

^C03  Ht  - 

♦ 

2L12  1 

^coa  Ot  -  £-\ 
v.  Y3/ 

.  arin^  A 1 1 

4- 

hi  “ 

a2  ( 

CBN 

F1DENTIAL 

sin2i'lT  i 

^~rT\ 

These  formulas  hold  for  any  0y{s,t). 

In  Ref.  1,  explicit  expressions  were  obtained  for  A^,  A22  for  the 

cone: 

-ft,  8*t 

(II. 7)  0v(s,t)  -  y2  +  p2e 


The  results  uore 
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12 


1  ~ 


*2  / 


/T  ♦  5  sin  2  .  /T 


ra  _ 

,  -7*-  ( 1  -  cos  I  L  T )  ♦  sin  ->  L  T  | 


2  sin  Sit  +  1  +  cos  i  .1  T I 


JA 


2  +  a,T  ♦  — - 


(II. Bill) 


9? 


1 


t  f1  +=j5, 

V  C2  / 


sin 


2  fit  +  1  +  COS  2  /It 


8in  / +  1  +  cos 


fit 


— 

2  *  V  +  77 


$ 


In  Figs.  1-7  are  shown  curves  of  the  quantity 


-  a2|  s-t| 


1/2 


&b  a 


'x(T)J 

Each  curve  is  determined 


function  of  T  for  07(s,t)  °  Yg  * 
by  specification  of  the  six  statistical  parameters  Y2>  Ty  ^ *  Pg* 

The  following  table  gives  the  parameter  values  associated  with  each 
curve £ 

One  word  of  explanation  about  the  table  is  in  order:  an  interesting 
case  is  the  case  in  which  the  velocity  dial  error  contains  a  "white  noise" 

P2  , 

component.  This  can  be  obtained  by  letting  a^-)  and  — - > /Cg  in  Eq.  (II. 7). 

The  exponential  component  of  the  velocity  dial  error  then  approaches  white 
ncise  with  spectra],  density  k  per  cycle. 
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T:  ’lie  following  tabic,  the  units  for  the  varic'is  parameters  are 
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ae  foil  ova: 


Y^  i  (nautical  miles)*1  (hours)  ^ 

-2 


t  (nautical  miles  )*  (hours) 

2 

Y-,  i  (nautical  miles’' 

2  **2 

Yj^  i  (nautical  miles )£  (hours) 
i  (hours)-'1 

^2  ’■  (nautical  miles )^  (hours) 

The  individual  curves  in  Figs.  1-7  ere  identified  by  circled  numbers • 
The  parameter  values  associated  with  each  curve  are  as  follows t 
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As  )  n  Case 
points  in  (0,T) 
Lot 

(III.1) 

Also  let 


(in.?) 


and 


(III. 3) 


The  quantities 


ill .  CASE  2:  ”  Y, 


,  let  U  -  tQ  •••<tn_:j  ”  T  be  n  equally  spaced  tine 


!_VI 


V 


-1 


(matrix  inverse) 


N, 


rM 


n-1 


/ _ 


[Vln 


iJ"o 


n-1 

[>]„  ’ia  ‘i  l) 

ij-o 


An  -  1110  M 

mw  L  Jn 

*12  ‘  ^  M 

n-*oo  L-  Jn 

n  lim  [Agj J 

n  — >00  L  jn 

'll'  A-i 2»  A22  31,0  functions  of  T. 
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:e!’eit»nce  1  gave  formulas  :.qs .  17. ly,  IV. 13,  H.17  of  Kef.  1)  fron 

TZ - 

whin.  could  be  derived  the  ferula  lor  OA(tjr).  The  rssult  can  be  expressed 
as  foil  Oils: 

Let 

*12*  hi 

(III. a)  y[ 

*22  “  hi  *  ~2  *  \ 

\ 

\  ^  °  hi  *22  ~  hi 

V.  > 

Then t  for  0  ■  t  *  T , 

_  r 

(in. 5)  c|( t,T)  -  j|  \2  -  2^12  t  ♦  t*> 

and  therefore 

r? —  if  _2 . 

(III.5-)  -(T)  *  ^|K22  ‘  2  K12  T*  R1I  ^  i 

Thefle  fonanlas  hold  for  any  0g(o,t). 

In  Ref.  1,  explicit  expansions  ware  obtained  for  hi*  h 2  lor 

the  case 

-ejs-tl 

(in. 6)  0p(3,t)  -  rx  * 
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’T  :  remits  -r.ia 


i 


An 


i'Eu(2-*a1T) 


1  - 


j 


(111.711) 


/}ll(2-a1T) 


T  I 


12 


~W{ 


1  - 


1  ♦ 


“ ^r" 


i,  ! 


(IIU7J.il) 


Ar 

~w 


1  + 


”lT.  1  fc4”!1' 

T  oTT 


1‘  1+  **T 

1  +  V7^T7 


In  Figs.  8-10  are  shown  curves  of 


c"2“" 

x(T) 


1/2 


as  a  function  of  T 


-uJs-t| 

for  0  (s.t)  “  y-  +  p.e  .  Each  curve  is  determined  by  specification 

B  1  rx 

of  the  six  statistical  parameters  Y2»  Yj(  Y^»  Op  Pp 

An  interesting  case  is  the  case  in  which  the  accolorometer  dial  error 

contains  a  "white  noise"  component.  This  can  be  obtained  by  letting 
Pi  v 

and  —  ■*  /t.  in  Eq.  (ill. 6) .  The  exponential  component  of  the  accelerometer 
®1  1 

dial  error  then  approaches  white  noise  with  spectral  density  per  cycle. 
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In  the 
follows ; 


■  tarile,  the  units  for  the  various  parametsrr  axe  as 


i  (nautical  stiles  )*"  (hours)  ^ 

2  _p 

Yj  i  (nautical  miles)  (hours ) 

o 

Yjj  «  (nautical  "'.■'Ire)’ 

Yjj  >  (nauti miles)1"  (hours  )“^ 

j  (hours)-1 

i  (nautical  miloa;”  (,  hoars)-11 


Tiie  pararoater  velues  associated  i,„th  each  curve  in  Figs.  8-10  are 


as  follows: 
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